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ABSTRACT 

We study the effect of a prolonged epoch of reionization on the angular power spec- 
trum of the Cosmic Microwave Background. Typically reionization studies assume a 
sudden phase transition, with the intergalactic gas moving from a fully neutral to a 
fully ionized state at a fixed redshift. Such models are at odds, however, with detailed 
investigations of reionization, which favor a more extended transition. We have mod- 
ified the code CMBFAST to allow the treatment of more realistic reionization histories 
and applied it to data obtained from numerical simulations of reionization. We show 
that the prompt reionization assumed by CMBFAST in its original form heavily con- 
taminates any constraint derived on the reionization redshift. We find, however, that 
prompt reionization models give a reasonable estimate of the epoch at which the mean 
cosmic ionization fraction was « 50%, and provide a very good measure of the overall 
Thomson optical depth. The overall differences in the temperature (polarization) an- 
gular power spectra between prompt and extended models with equal optical depths 
are less than 1% (10%). 



1 INTRODUCTION 



At a redshift of z ~ 1100 the intergalactic medium (IGM) 
recombined and remained neutral until the first sources of 
ionizing radiation formed. While the distribution and evo- 
lution of these sources are unknown, the overall process of 
IGM reionization is fairly well understood, and can be di- 
vided into three phases. First, individual HII regions devel- 
oped around the sources. Then, these ionized regions grew 
in number and size until they overlapped, producing a sud- 
den increase of the photon mean free path. Finally, when all 
underdense regions and voids were completely ionized, pho- 
tons penetrated into overdense clumps and filaments, bring- 
ing the reionization process to completion at a "reionization 
redshift." 

It is clear that the mystery of the nature and evolution 
of the ionizing sources is intimately tied with the time scale 
over which these phases took place. Thus there is a great 
deal of physical information that could be gathered if the 
reionization redshift, Zi, and its duration, Az, can be firmly 
established. 

Apart from the study of the (HI and Hell) Gunn- 
Peterson effect, which has not yet yielded conclusive re- 
sults (Becker et al. 2001; Gnedin 2001; Theuns et al. 2001), 
mcEisurements of Cosmic Microwave Background (CMB) 
anisotropics have the greatest potential for constraining 
these important quantities (Griffiths et al. 1999). The scat- 
tering of CMB photons by free electrons damps the angu- 
lar power spectrum of primary anisotropics by a factor of 
e~^^ for large angular multipoles I ^ 100 (Tegmark & Zal- 
darriaga 2000), where r is the Thomson optical depth. Us- 



ing the currently available data, consistency with the lack 
of Gunn-Peterson trough and with the observed peak in 
the angular power spectrum at ^ ~ 200 (De Bernardis et 
al. 2000; Hanany et al. 2000; Padin et al. 2001) is able to 
constrain 0.02 < r < 0.44 (De Bernardis et al. 1997; Grif- 
fiths et al. 1999; Tegmark et al. 2001; Griffiths & Liddle 
2001), although these results are somewhat dependent on 
the cosmological model assumed. Much better constraints 
are expected from polarization studies to be carried out by 
future satellites such as S'POriQ and PLANCt^. 

A second method of quantifying reionization is to ex- 
amine Zi, and several recent studies have attempted to use 
the CMB to derive this redshift directly, using codes that as- 
sume a sudden epoch of reionization (e.g, Hu et al. 95; Seljak 
& Zaldarriaga 1996). For example, Schmalzing et al. (2000) 
use MAXIMA data in combination with cosmological pa- 
rameters from independent measurements of Big Bang Nu- 
cleosynthesis and X-ray cluster data to constrain Zi. By per- 
forming a analysis, they conclude that Zi > 15(8) at the 
68% (95%) confidence level. Similarly, Naselsky et al. (2001), 
in addition to providing support to Schmalzing et al. results, 
showed that polarization spectra are very sensitive to the 
reionization redshift. 

An important caveat is implicit in these applications 
however, namely the assumption of prompt reionization. 
While taking Az = is a reasonable first step, more de- 
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tailed simulations show that this approach paints a picture 
of reionization in only the broadest of strokes. This is par- 
ticularly true as the onset of reionization raises the overall 
temperature of the IGM, suppressing the further formation 
of objects that are too weakly bound gravitationally to over- 
come the drastic increase in thermal pressure (Barkana & 
Loeb 1999). Thus recent numerical simulations of reioniza- 
tion (Ciardi et al. 2000, hereafter CFGJ; Bruscoh et al. 2000; 
Gnedin 2000; Benson et al. 2000) have shown that the evolu- 
tion of the mean ionization fraction is slow, has a nonlinear 
dependence on redshift, and occurs in an extremely patchy 
manner. All of these details may considerably affect the de- 
termination of Zi. 

In this paper, we aim to clarify the effects of realistic 
reionization scenarios on the observables typically used to 
examine the reionization epoch. To this end we have modi- 
fied the most heavily relied on theoretical code for comput- 
ing CMB fluctuations, CMBFAST (Seljak & Zaldarriaga 1996), 
to allow the treatment of such reionization histories, helping 
to clarify the best approach in trying to quantify and define 
the epoch of reionization. 

The structure of this work is as follows. In §2 we dis- 
cuss numerical models of reionization and how these have 
been incorporated into the CMBFAST code. In §3 we apply 
these techniques to study the impact of these models on the 
temperature and polarization spectrum of the CMB, and 
conclusions are given in §4. 



Xe{z) are obtained from high-resolution cosmological simula- 
tions within a periodic box of comoving length L = 2.55/i~^ 
Mpc. The source properties are calculated taking into ac- 
count a self-consistent treatment of both radiative feedback 
from ionizing and H2-photodissociating photons and stel- 
lar feedback regulated by supernovae from massive stars. 
There are two main free parameters in the simulations: the 
fraction of total baryons converted into stars, fbi,, and the 
escape fraction of ionizing photons from a given galaxy, 
/esc- A critical discussion of these parameters is given in 
CFGJ, and here we examine four different models, fixing 
(/b* = 0.012, /esc = 0.2) in run A, (0.004, 0.2) in run B, 
(0.15, 0.2) in run C , and (0.012, 0.1) in run D. As run A, in 
which Zi ^ 11, gives the best agreement between the derived 
evolution of the cosmic star formation rate and observations 
at z ^ 4 (Steidel et al. 1998, CFGJ), we choose this as our 
fiducial model, and rely on runs B, C, and D to study the 
effects of model uncertainties on our conclusions. 



3 PROMPT VS. SIMULATED REIONIZATION 

Equipped with this more general tool, we now quantitatively 
contrast the prompt reionization Xe,p{z) with the one de- 
rived from the simulations, Xe,s[z). As CMBFAST can be run 
either by assigning Zi or r, we are able to examine both these 
cases separately. 



2 REIONIZATION HISTORY FROM 
SIMULATIONS 

To assess the effects of more physical reionization histories 
on the determination of Zi, we have modified the Boltzmann 
code CMBFAST such that the redshift evolution of the volume- 
averaged mean hydrogen ionization fraction, Xf,[z), can be 
specified, reading in data from more detailed simulations of 
reionization.^ Our modified version of CMBFAST is then able 
to consider three types of reionization histories: (i) prompt 
reionization models in which the total optical depth, r, is 
chosen (n) prompt models in which the reionization redshift, 
Zi, is fixed, and [Hi) models in which x^{z) is read from an 
external file. 

As in Bruscoli et al. (2000), we consider a set of reion- 
ization simulations studied in CFGJ, which model reioniza- 
tion from stellar sources, including Population III objects. 
These simulations were developed in a critical density CDM 
universe with cosmological parameters h = 0.5, ag = 0.6, 
and Q.b = 0.06 (where h is the Hubble constant in units of 
100 kms~^ Mpc~^, (Tg is the rms variation of density pertur- 
bations at the 8/i~^ comoving Mpc length scale, and 0,^, is 
the baryonic density in units of the critical density) but span 
a large range of model parameters, and thus are suitable for 
drawing conclusions as to the effect of extended reionization 
in more general cosmologies. 

The spatial distribution of the sources and ionized re- 
gions at various cosmic epochs and the overall evolution of 



■t A copy of this version of CMBFAST is publicly avai lable at 
tittp : //wwH. arcetri . astro . it/scieiice/cosmology/CMB| 



3.1 Assigning the Reionization Redshift 

In Fig. we plot the redshift evolution of the hydrogen 
ionization fraction and the Thomson optical depth for the 
prompt and simulated cases, choosing in the prompt case 
the same value of Zi — 10.9 as deduced from run A. In this 
plot, the smoother, nonlinear increase of Xe,s{z) is evident, 
with Az ~ 15. Note also that the optical depth for the two 
reionization histories are quite different, and that Ts — 0.08 
for run A and Tp — 0.04 for prompt reionization. If we define 
the discrepancy between the prompt and simulated visibility 
functions g{z) as 

Afl _ \gs{z) - gp{z)\ 

where g{z) is the probability that a photon reaching the 
observer last scattered in the redshift range z and z + dz, 
then we find that /S.g/g « 5% for a large range of redshifts 
between recombination and the beginning of reionization 
(25 < 2 < 1000). 

The most important comparison between the prompt 
and fiducial cases, however, is made in terms of the dif- 
ferent resulting CMB temperature and polarization angular 
power spectra, Ci and Pi, as these are the observable quan- 
tities. These are shown in the upper left panel of Fig. 2, in 
which differences are particularly evident in the polarization 
spectra, which are more sensitive to the reionization history. 
Reionization introduces a characteristic bump in the lower 
multipoles (^ = 5 — 10) that tends to shift to higher ^'s when 
r is increased. The amplitude of the bump in the Pi spec- 
trum decreases with t and therefore the signal is weaker for 
the prompt reionization case. In particular the bump am- 
plitude in the angular spectrum Pi is roughly proportional 
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Figure 1. H ionization fraction, Xe, and Thomson optical depth, 
T, versus redshift for run A (solid lines) and for prompt reioniza- 
tion (dotted lines) when CMBFAST is run choosing Zi = 10.9. The 
dots represent the output from the simulation. 

to (for r ^ 1), while its position scales as £peak oc 
(Zaldarriaga 1997; Fabbri 1999). 

To highlight these differences, we plot the normalized 
discrepancy of temperature and polarization angular spec- 
tra [defined analogously to eq. (1)] between the models in 
the lower left panel. Here find that ACe/Ce < 0.07 while, 
around I ~ 10, APi/Pi ~ 1. These results make it clear 
that although the prompt and simulated reionization his- 
tories share the same value of Zi, they produce noticeably 
different angular power spectra, and thus correct constraints 
on the reionization redshift cannot be obtained by equating 
Zi in the prompt and simulated runs. 




10 102 IQ3 10 102 IQ3 
1 1 

Figure 2. Top left: temperature (Ci) and polarization (Pi) an- 
gular power spectra versus the multipole £ for run A (solid line) 
and for prompt reionization (dotted line) when CMBFAST is run 
with Zi = 10.9. Bottom left: the normalized temperature and po- 
larization discrepancy as a function of angular multipole £ for this 
case. Here the solid and dotted lines refer to Ci and Pg respec- 
tively. Top right: temperature (Ci) and polarization (Pi) angular 
power spectra versus the multipole £ for run A (solid line) and 
for prompt reionization (dotted line) when CMBFAST is run with 
r = 0.08. Bottom right: the normalized temperature and polariza- 
tion discrepancy as a function of angular multipole £ for this case. 
Here the solid and dotted lines refer to C( and Pg respectively. 



3.2 Assigning the Optical Depth 

An alternative approach is to compare our fiducial run with 
a prompt model that differs in its Zi but has the same Thom- 
son optical depth, r — 0.08. This results in a value for the 
prompt reionization redshift of Zi^p = 17.1. Note that this 
value is much larger than the simulated redshift of reioniza- 
tion Zi^s = 10.9, and corresponds to the redshift at which 
the simulated ionization fraction is only Xe,3 ~ 0.5. 

Comparing these two runs we find, however, that the 
discrepancy between the two visibility functions is only 
Ag/g ~ 10~*. This small change in the visibility function 
results in an equally small change in the observed CMB 
anisotropics. This can be seen on the top right panel of in 
Fig. ^ in which the angular power spectra are almost in- 
distinguishable. Indeed, the fractional discrepancy between 
these runs is only ACi/Ce < 0.01 and APijPi < 0.06 as 
shown in the bottom right panel in this figure. 

Thus the two different reionization models produce 
comparable angular power spectra if r is specified rather 
than Zi . In fact if we try to recover the reionization redshift 
from the standard formula Zi = 8.9{T/hQ,i,)^^''Qg''^ (Peebles 
1993, Tegmark et al. 2000) we get Zi = 17.1, again a redshift 
at which a;e ~ 0.5. It is clear then that by analyzing CMB 



power spectra using a prompt model, one can draw reason- 
able conclusions as to the overall Thomson optical depth, 
whereas the reionization redshift is much more uncertain. 

3.3 Quantification of Errors & Tests of Our 
Approach 

In order to quantify the difference between the fiducial 
model and the prompt model with the same optical depth 
further, we plot in Fig. ^ the quantity 

/(£^„.)=5^^(^)^f... (2) 

2 

as a function of the multipole number £max assuming full 
sky coverage, fsky = 1- For multipole values greater than 
100, exceeds unity, the value corresponding to the cos- 
mic variance error of 2C| / {2£ + 1). This means that one can 
find a statistically significant difference between the prompt 
and fiducial models only by looking at angular multipoles 
with £ ^ 100. However, in Fig. ^ we saw that the largest dif- 
ferences between the two histories occur at lower multipole 
values (i* = 20 — 40). Thus although these reionization his- 
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Figure 3. of C** (solid lino) and of Pi (dotted line) versus 
^max (see eq. 2). 



tories are considerably different, it is almost impossible to 
discriminate between them observationally. 

Finally, we extend our analysis to the runs B, C, and 
D of CFG J, comparing them with prompt runs with equal 
values of Ts- These results are summarized in Table 1. As in 
run A, the values of Zi^p found for these models correspond 
to times at which le.s ~ 0.5, when the IGM was in the midst 
of changing from a neutral state to an ionized one. 

Is is clear from this table that it is not possible to 
put a constraint on the reionization redshift directly from 
CMBFAST. In fact, as can be easily seen in Table 1, the differ- 
ences between Zi^s and Zi^p cannot be easily parameterized 
because they depend on several effects that influence the du- 
ration Azs of the reionization process. These effects include 
radiative and stellar feedbacks, the star formation efflciency, 
the photon escape fraction, the spatial and mass distribu- 
tion of the ionizing sources, and even the dumpiness of IGM, 
which is beyond the scope of the CFGJ simulations. 

We have conducted a number of checks and convergence 
tests to assess the robustness of our approach. Our first check 
was to examine the overall, numerical value of the optical 
depth, which differs slightly from the input value due to the 
finite integration time step. Here we found that the output 
TS differ by less than 5 x 10~^ from the specified values. 
This corresponds to ACe/Ci Si 10~*, i.e. 30 times smaller 
than the actual difference between the simulated and prompt 
runs. 

Next, we carried out a series of convergence tests, to 
ascertain the effects of fc-mode sampling on our results. 
Here we restrict our tests to run A for which we increased 
the k resolution by 50%. A comparison between the origi- 
nal and resampled simulated runs showed that at most an- 
gular scales ACe/Ce and APe/Pe are within 0.1%; larger 
differences (up to ~ 0.7%) are found at higher multipole 
numbers {£ ^ 700). Comparing the resampled run A with a 
resampled prompt run with equal input optical depth, we 
found ACe/Ce and APi/Pe values that were virtually in- 
distinguishable from those obtained with low k resolution 



Finally, we examined the effect of time sampling at both 
k resolutions, comparing two prompt runs with equal optical 
depth and k resolution, but forcing one of them to the same 
time step as in simulated run A. Again, we found ACe/Ci 
and APt/Pe ^ 10~*, far smaller than the discrepancies be- 
tween the prompt and simulated runs. 



4 CONCLUSIONS 

One of the first stages of nonlinear structure formation, 
reionization marked an important transition from a dark 
and relatively simple universe to one filled with a dazzling 
array of stars, galaxies, quasars, and other nonlinear ob- 
jects. And although one of our best probes of this transition 
is through the measurement of CMB fiuctuations, the pro- 
cess of reionization itself is much more dependent on the 
complicated astrophysical issues important at low redshifts 
than the linear issues important at 2 ~ 1100. 

In this work, we have explored this transition, quan- 
tifying the impact of realistic scenarios of reionization on 
the angular power spectrum of the Comic Microwave Back- 
ground. While standard estimates assume prompt reioniza- 
tion, we have considered instead a range of simulated mod- 
els, each with a prolonged reionization epoch. We find that 
equating the redshift of full IGM ionization between these 
simulations and models that assume instantaneous reioniza- 
tion leads to widely discrepant temperature and polarization 
spectra. On the other hand, equating prompt and extended 
models with the same overall optical depth leads to differ- 
ences in anisotropics that are nearly undetectable. In this 
case the redshift of complete ionization is lost in the compli- 
cated details of the phase transition, and comparisons yield 
Zi values corresponding to roughly the point of 50% ioniza- 
tion in the simulations, although even this value is model de- 
pendent. It is clear then that while Zi is useful as schematic 
tool, it is the total optical depth that is most accurate in 
providing a definition of the reionization epoch. 

MB thanks M. Zaldarriaga for discussions and hospi- 
tality at IAS, Princeton where part of this work has been 
carried out. ES has been supported in part by an NSF MPS- 
DRF fellowship. 
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Run 




z- C") 




2- 








A 


15 


10.9 


0.080 


17.1 


0.46 


^ 0.01 


;S 0.06 


B 


17 


8.3 


0.059 


13.9 


0.51 


< 0.007 


< 0.1 


C 


15 


15 


0.098 


19.8 


0.48 


^ 0.01 


;$o.i 


D 


17 


8.3 


0.064 


14.7 


0.47 


;$ 0.008 


;$o.i 



Table 1. ("'Duration of reionization; ('''Redshift of complete 
reionization (simulated); ('^'Thomson Optical Depth (simulated); 
('''Redshift of complete reionization (prompt); t*^' Hydrogen ion- 
ization fraction at p (simulated). 
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